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ABSTRACT 

By comparing sensitive Spitzer Space Telescope (SST) infrared and Green Bank Telescope 21 cm 
observations, we are able to report the first detection of dust emission in Complex C, the largest High 
Velocity Cloud in the sky. Dust in the region of Complex C studied here has a colder temperature 
(T = 10.7io'g K - la) than the local interstellar medium (T = 17.5 K), in accordance with its 
great distance from the Galactic plane. Based on the metallicity measurements and assuming diffuse 
Galactic interstellar medium dust properties and dust/metals ratio, this detection could imply gas 
column densities more than 5 times higher than observed in Hi. We suggest that the dust emission 
detected here comes from small molecular clumps, spatially correlated with the Hi but with a low 
surface filling factor. Our findings imply that the HVCs' mass would be much larger than inferred 
from Hi observations and that most of the gas falling on the Milky Way would be in cold and dense 
clumps rather than in a diffuse phase. 
Subject headings: Galaxy: Halo, ISM: Clouds, ISM: Dust 



1. INTRODUCTION 



Sin ce their discovery in HI observations l)Muller et alJ 
1196^1 . High Velocity Clouds (HVCs) have been the tar- 
get of numerous studies (see Ivan Woerden et all (|2004T) 
for a thorough review) but nevertheless remain puzzling. 
It is now widely considered that many of them might 
be infalling clouds fueling the Galaxy with low metallic- 
ity gas. This hypothesis received strong support from 
ultraviolet observ ations showing that H VCs have a sub- 
solar metallicity jWakker et alJ 11999^1 and a D /H ra- 
tio compatible with p rimordial abundances of deuterium 
l|Sembach et alJl2004jl . 

The present study focuses on the Spitzer Extra- 
Galactic First Look Survey (XFLS) field, a diffuse HI 
area (Nhi ~ 2 x 10 20 cm" 2 ) at high Galactic latitude 
(/ = 88°. 3, b = 34°. 9) programmed to st udy galaxy popu- 
lations and clustering ijFang et alJ2004|) . The XFLS field 
is located on the edge of Complex C, a high- velocity gas 
structure spanning more than 1500 square degrees on the 
sky. Complex C h as a subsolar metallicity (~ 0.1 — 0.3) 
l|Tripp et all 12003). it contains ionised gas observed in 
Ha emission ijTufte et aD [19981 and O VI absorption 
l)Sembach et al.l2003t ) and i s located at a d istance greater 
than 5 kpc from the sun Dust emis- 

sion has be en unsuccesfully looked for in HVCs using 
IRAS data l|Wakker & Boulangerlfl 986ft . Upper limits 
were found to be compatible with their low Hi column 
densities, low metallicity (and therefore low dust/gas ra- 
tio) and their large distance (i.e. fainter radiation field). 
In this study we revisit the dust/gas correlation in HVCs 
using new infrared and Hi data. 

2. DATA 
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In this study we used Spitzer Space Telescope (SST) 
MIPS 24 and 160 data of the XFLS. The field 
was covered by MIPS in 8 scan maps that were com- 
bined and destriped, then regridded to 16" pixels and 
median-smoothed. We complemented the SST data 
with improved IRAS 60 and 100 /im maps (IRIS - 
iMiville-Deschenes fc Lagachel ©05)). The FWHM an- 
gular resolution of the IRIS maps is ~ 4' with a pixel size 
of 1'.5. The calibration uncertainties on diffuse emission 
of the SST and IRIS maps is 10 %. Bright point sources 
were removed from all infrared maps (SST and IRIS) 
before convolution by a Gaussian beam to bring all of 
them to the Green Bank Telescope (GBT) 9'. 2 FWHM 
resolution. Fig. Q] shows the convolved maps. 

The Hi (21 cm) data we used in this study were ob- 
tained with the G BT, previously published an d kindly 
made available by lLockman fc Condon! l)2005l) . These 
observations cover a 3° x 3° area centered at (J2000) 
a = 17' l 18 m , S = +59°30', which encloses the SST 
observations. The GBT data were corrected for stray- 
radiation contamination. The sampling was done at 1'5, 
the FWHM beam is 9.2', the channel width is dv =0.51 
km s" 1 and the noise level in a single channel is 0.05 K. 

3. HI DECOMPOSITION 

Fig-Elpresents the average 21 cm spectrum of the field 
and two typical lines-of-sight. In addition to Complex 
C, the XFLS 21 cm data revealed two Intermediate- 
Velocity Clouds (IVC1 and IVC2) and a local compo- 
nent typical of high-latitude Hi emission with a narrow 
cold co mponent superimposed on a broad warm com- 
ponent l|Lockman fc Condon! 12005ft . In this section we 
estimate the column density of each of these Hi compo- 
nents. The HVC is well isolated in velocity space and 
its integrated emission was computed by adding all the 
channels where —210 < v < —132 km s _1 . On the other 
hand the two IVCs and the local component overlap in 
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Fig. 1. — Infrared maps of the XFLS field at 24 and 160 micron 
from Spitzer and 60 and 100 micron from IRIS. All maps were 
convolved by a Gaussian beam to bring them to the GBT 9.2' 
FWHM resolution. 
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Fig. 3. — Hi column density of the four Hi components (Local, 
HVC, IVC1 and IVC2). All maps are in units of 10 20 cm" 2 . An 
opacity correction was applied to convert the 21 cm integrated 
emission to Hi column density - see text for details. 
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Fig. 2. — Decomposition of the 21 cm data. The bottom curve 
shows the 21 cm spectrum averaged over the whole field. The 
middle and top curves (shifted for clarity) are representative 21 
cm spectra where the IVC1 (top) and IVC2 (middle) components 
are observed. The shaded rectangles indicate the velocity range 
used to extract the Hi components. 



velocity. The integrated emission of the two IVC compo- 
nents were computed by adding all the channels between 
-40 < v < -14 km s" 1 for IVC1 and -50 < v < -29 km 
s _1 for IVC2 from which we subtracted the local compo- 
nent contribution estimated by interpolating a baseline 
on each spectrum (see Fig.|5J). The local component's in- 
tegrated emission was computed by adding all channels 
between —48 < v < +25 km s _1 from which the two 
IVC contributions were removed. There is still contri- 
bution from the local components in a few channels at 
v < —48 km s _1 but they were not added, as the emission 
is very diffuse and its structure within the noise. 
In the optically-thin case (i.e. when the 21 cm bright- 



ness temperature is much smaller than the Hi kinetic 
temperature), the 21 cm integrated emission is propor- 
tional to the Hi column density: Nm(x,y) — 1.823 x 
10 18 J2 V Tb{x, y, v)dv. The local component is domi- 
nated by a broad (i.e. warm and optically-thin) compo- 
nent but a fraction of the emission comes from cold gas 
(narrow features) that needs to be corrected for opacity 
to properly estimate the column density. To do so we es- 
timated the fraction of warm and cold gas on each line-of- 
sight using a Gaussian decomposition of each spectrum. 
For each Gaussian component we estimated its opac- 
ity correction, considering that thermal and turbulence 
broadening equally contribute to the velocity dispersion 
a v (i.e. spin temperature is T s = 121g, 2 72), instead of 
the u sual isothermal assumption ijLockman fc Condonl 
( 2005) used T s = 80 K). This allowed us to compute 
a larger correction for the narrow than for the broad 
components. Overall the opacity correction for the local 
component is lower than 6 %. It reaches its maximum 
value in the bright filament (south-western part of the 
field, see Fig. |3J where our estimate of the Hi kinetic 
temperature gets as low as 50 K. The opacity correction 
for the two IVCs and the HVC is less than 1% in ac- 
cordance with their broad lines ( FWHM> 8 km s _1 ). 
Fig. |3 presents the Hi column density maps of the four 
Hi components. 

4. INFRARED-HI CORRELATION 

The main goal of this study was to separate and deter- 
mine the infrared emission of the HI components identi- 
fied with the 21 cm observations. To do so we considered 
that the infrared map I\(x,y) at wavelength A can be 
represented by the following model 



h{x,y) =^a\Nj II (x,y) + C x {x,y) 



(1) 



where N l HI (x, y) is the column density of the i HI com- 
ponent, a\ is the infrared-Hi correlation coefficient of 
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component i at wavelength A and C\{x,y) is a residual 
term. 

To estimate the correlation coefficients a\ we per- 
formed a chi-square minimisation on the Hi and infrared 
data based on the model given by Eq. ^ The infrared- 
Hi correlation coefficients are collected in Table ^ an d 
shown in Fig. 0] The standard deviation of the residual 
maps C\(x,y) are also given in Tabled 

The uncertainties on a\ given in Table ^ include cali- 
bration and statistical uncertainties, taking into account 
the data sampling. To estimate the impact of our opacity 
correction we made the same analysis with optically-thin 
Hi maps and obtained the same results within the un- 
certainties. 

The infrared-Hi correlation presented here shows that 
the model used (Eq. ^| for the infrared emission is 
satisfying. At each wavelength the standard deviation 
of the residual ter m C\(x,y) is in go od concordance 
with the expected llLagache et al.l 120031) and measured 
l)Miville-Deschenes et al.ll2002fl levels of the Cosmic In- 
frared Background (CIB) fluctuations (see Tabled, ex- 
cept at 160 /im where we report a significant residual 
excess. As the residues at 60 and 100 are in ac- 
cordance with the CIB level, this analysis also leads us 
to conclude that dust associated with the Warm Ionised 
Medium (WIM) l|Lagache et al.lll999]) does not signifi- 
cantly bias our decomposition of the infrared emission in 
this field. 

The striking result of this analysis is certainly the 
significant detection of the infrared emission associ- 
ated with an HVC at all wavelengths. Our result at 
100 um. is compatible wit h the upper limit given by 
lLockman fc Condon! Ipffll . 

5. DUST PROPERTIES OF COMPLEX C 

Using the infrared-Hi correlation coefficients we now 
want to put contrains on the properties of the dust, es- 
pecially in the HVC. 

We assume that the emission from big grains at ther- 
mal equilibrium with the radiation field is a modified 
black body : 

h = t x B u (T bg ) (2) 

where £?„ is the Planck function, Tbg is the big grains' 
equilibrium temperature and t\ the opacity, often ex- 
pressed as t\ = Nhi^x where t\ is the dust emissivity 
per H atom. The emissivity is generally described as 
a power law : e\ = £qA - ^ where (3 is the emissivity 
spectr al index. Based on FIRAS data, iBoulanger et al.l 
( 1996) give the following expression for the big grains' 
emissivity law in the local interstellar medium: 

(3) 

In this scheme, and at wavelengths where the emis- 
sion is dominated by big grains, the infrared-Hi correla- 
tion coefficients can be expressed as a\ = €\B„(Tbg)- 
Assuming that Very Small Grains (VSGs) do not con- 
tribute sig nificantly i n the far-infrared and taking (3 = 2 
(praine & Leelll984l IBoulanger et alJll996j) . we use the 
100 and 160 /im infrared-Hi correlation coefficients to es- 
timate the big grains' temperature T BG of each Hi com- 
ponent 

c4 _ B,(10Q^m,T BG )16Q 2 
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Fig. 4. — Hi -infrared correlation coefficients for the Local, 
HVC, IVC1 and IVC2 components shown in Fig. |3] See Table [I] 
for details. 
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and the dust emissivities at 160 /im : e{ 
a\ 60 / B„(160[j,m, T l BG ). The big grain temperature and 
the 160 /im emissivity of each Hi component are given 
in Table |3 In this table are also given the corresponding 
1 and 3<T ranges on T BG and e\ 60 computed using 1 and 
3a uncertainties on a\ 00 and a\ 60 . 



6. DISCUSSION 



The values of T l ^ al 



17.4±^ and = 3.0±^ x 

10~ 25 cm 2 per H atom found for the local component 
are in very good agreement (within la uncertainties) 
with the average values f o und fo r the Solar neighbor- 
hood by Boulang er et al.l l)1996|) using COBE data : 
T B ° G n = 17.5 K and efg z n = 2.4 x 10" 25 cm 2 per H 
atom. This is an important validation of the photomet- 
ric accuracy of the SST-MIPS data. On the other hand 
dust in the HVC is found to be significantly colder than 
in the local interstellar medium, even at the 3a limit 
{T BG C = 10.71 2 q), which is consistant with a lower ra- 
diation field than in th e Solar neighborhood due t o its 
distance to the Galaxy l)Wakker fc Boulangerlll98ri) . In 
addition we report an emissivity eigo higher in the HVC 
than in the local ISM (see Table [2J. An eventual contri- 
bution of VSGs to the far-infrared emission would have 
the effect of lowering Tbg an d increasing ei6o even more. 

For dust properties typical of the local interstellar 
medium, e\ should be proportional to the dust-to- 
gas mass ratio which, for standard metal depletion on 
grains, would scale with the metallicity. The fact that 

C HVC / c soln 
^60 



/ e i60 > fia) is greater than the metallic- 



n60 



5^(160^, T|; G )100 2 



ity in Complex C (Zhvc/Zq = 0.2 ± 0.1 according to 
Trip D*et alJ | 2003^ could indicate a dust (and therefore 
gas) column density in Complex C more than 5 times 
larger than given by the Hi emission. 

Could it be dust associated with ionised gas ? Ha 
and other ionised atomic lines have been detected in sev- 
eral HVCs. It is expected to come from a photoionized 
layer at the surface of the clouds which plunge into the 
Galactic halo. The column densit y of such a layer is 
2 - 4 x 10 19 cm- 2 in Complex C (IWakker et all 1999) 
which is similar to the Hi column density in the HVC 



4 



TABLE 1 
Infrared-Hi correlation 



A 


"a 


a Hva 


a ivoi 


a{ vo. 


a (Residue) 


CIB estimate 


H 










(MJy sr" 1 ) 


(MJy sr" 1 ) 


24 


0.017±0.002 


0.007±0.002 


0.039±0.006 


0.014±0.005 


0.004 


0.004 


60 


0.16±0.02 


0.05±0.01 


0.35±0.04 


0.31±0.04 


0.03 


0.02 


100 


0.80±0.08 


0.055±0.015 


0.70±0.09 


0.9±0.1 


0.09 


0.07 


160 


1.7±0.2 


0.8±0.1 


2.7±0.4 


1.4±0.4 


0.30 


0.11 



Infrared-Hi correlation coefficients (ct^) given in MJy sr 1 (10 20 H atoms) 1 cm 2 for the four Hi components (Local, HVC, IVC1 and 
IVC2) at 24, 60, 100 and 160 fim. The uncertainties on each is the la uncertainty taking into account the statistical variance and 
instrumental uncertainties. Column 6 gives the standa rd deviation of the re sidual map and Column 7 is an estimate of the Cosmic Infrared 
Background (CIB) fluctuations based on the model of Lagachc ct al. 1 2003). 



TABLE 2 
Dust parameters 





Tbg 


la 


3a 


ei60 


la 


3a 


Local 


17.4 


16.3-18.8 


14.3-22.6 


3.0 


1.9-4.9 


0.60-12.5 


HVC 


10.7 


9.9-11.6 


7.7-13.6 


36.0 


16.8-83.7 


3.8-1442.8 


IVC1 


14.6 


13.6-15.8 


11.8-19.4 


13.2 


7.0-24.1 


1.6-79.6 


IVC2 


19.4 


17.1-23.1 


14.1-107.4 


1.5 


0.5-3.6 


0.003-15.5 



Big Grain temperature (TgQ, given in Kelvin) and dust emissivity 
at 160 fim (ei60 = tiqq/Nhi given in 10 — 25 cm 2 per H atom) and 
their corresponding 1 and 3 a ranges for each Hi component. The 
value s for the local interstel lar medium deduced from the COBE 
data (Boulanger et al. 1996) are T b q = 17.5 K and eieo = 2.4 X 
10 -25 cm 2 . 



seen in XFLS and therefore not enough to explain the 
high infrared emissivity measured here. In addition only 
the ionised gas spatially correlated with the Hi would 
contribute to the infrared emissivity. Could it be dust 
associated with diffuse H2 ? Two detections of H2 in ab- 
sorption with FUSE have been reported b ut the column 
densit ies are much lower than those of Hi ijRichter et al.l 
l2001a|L Other attemps have been unsucces ful; CO was 
not found in emission tfflakkexj^LaU llM ZD nor HCO + 
in absorption ijCombes & Charmandarisl l200C0 . Further- 
more abundance measurements in HVCs show no sign 
of iron depletion <|Murphv et al.ll2f)0fl iTrTpD et alJl2003) 
which has been used as a strong argument in favor of 
very low dust in HVCs. 

To reconcile our findings with previous work we sug- 
gest that HVCs could be multiphase gas including small 
and dense molecular clumps with a low surface filling 
factor (and therefore undetected with pencil beam ab- 
sorption observations) bathing in a diffuse phase seen 
in absorption. The dust emission reported here would 
mainly originate in these clumps, spatially correlated 
with Hi, which would have a larger iron depletion than 



in the diffuse phase of HVCs. Evidence of the exis- 
tence o f such small molec ular clumps was found in cirrus 
clouds | Hei1;haus er3l2f)04l) and in the Magellanic Stream 
l)Sembach et al. 120011) . Additional observations of HVCs 
with the Spitzer Space Telescope are needed to confirm 
our result. 

If our result reveals a general property of HVCs, it 
would mean that their mass is much larger than what 
can be inferred from Hi observations and that most of the 
gas is in a dense phase. An infall of 1 M^/year is needed 
to explain the star formation rate, the distribution of 
stellar metallicities with age and the presence of a bar in 
the Milky Way and in other galaxies. The current idea 
is that this infalling gas is diffuse (ionised and neutral) 
resulting from the cooling of very hot halo gas. Our 
result changes this picture by showing that most of the 
infalling gas is in dense and cold clumps, leaving open the 
question of the origin and stability of such structures. 

HVCs typically contribute to < 10% of the Hi column 
density at high latitude, but with its temperature lower 
than in the local interstellar medium, dust in HVCs could 
contribute a much larger fraction to the sub-millimeter 
and millimeter emission. Observations in this wavelength 
range, with Herschel and Planck for instance, could be 
used as a tracer of the column density of HVCs and as 
a distance indicator, using the dust temperature to put 
constraints on the radiation field strength. If the detec- 
tion presented in this study is confirmed and is typical of 
HVCs, it could also have a strong impact on the analysis 
of observational cosmology data (e.g. WMAP, Planck) 
that do not to this day take this foreground component 
into account. 

The authors thank the anonymous referee for a careful 
reading of the manuscript and very helpful comments. 
M.A.M.D. acknowledges support from The Canadian 
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Astrophysics Data System Bibliographic Services. 
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